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An  analytic  expression  is  derived  for  the  mean  cosine  of  the  Fournier-Forand1 2  phase  function.  This 
expression  and  the  power  law-  index  of  refraction  relationship  of  Mobley3  are  used  to  parameterize  the 
Fournier-Forand  phase  function  by  its  mean  cosine  in  a  similar  manner  to  the  Flenyey-Greenstein4 
function. 

Radiative  transfer  theory  makes  extensive  use  of  the  mean  cosine  of  the  phase  function  to  compute  the 
evolution  of  the  light  field  under  scattering.  In  this  paper  an  analytic  expression  for  the  mean-cosine  of  the 
Fournier-Forand  phase  function  in  terms  of  Hypergeometric  functions  is  derived.  This  expression  reduces 
to  a  simple,  easy  to  evaluate,  very  strongly  convergent  series  that  is  a  function  of  both  mean  index  of 
refraction  and  the  Junge  size  distribution  inverse  power  law  exponent.  By  judicious  use  of  the  relationship 
between  the  index  of  refraction  and  the  power  law  proposed  by  Mobley  it  becomes  possible  to 
parameterize  completely  the  Fournier-Forand  phase  function  in  terms  of  its  mean  cosine.  This  result  is 
more  complex  but  still  analogous  to  the  situation  that  prevails  with  the  I  Icnyc-Grccnstcitr  phase  function 
where  the  asymmetry  factor  is  identical  to  the  mean  cosine.  The  implications  of  using  this  expression  for 
the  mean-cosine  in  the  multiple  scattering  regime  of  radiative  transfer  are  explored  along  the  lines  recently 
suggested  by  Piskozub  and  McKee6. 

In  order  to  obtain  an  expression  for  the  mean  cosine  we  need  to  evaluate  the  following  two  integrals  in 
sequence: 
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We  therefore  need  to  first  derive  p[x,6 )  from  first  principles.  We  start  with  the  single  particle  phase 
function  approximation.  We  want  this  function  to  be  normalized  to  unity. 
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Note  also  that,  as  in  the  case  of  the  derivation  of  the  full  Foumier-Forand  phase  function,  we  approximate 
the  single  particle  Airy  function  angular  diffraction  pattern  by  the  following  approximation: 
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So  the  single  particle  normalized  phase  function  approximation  is 


f{x,6)=^~ 
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To  get  the  full  phase  function  we  then  need  to  integrate  over  all  x  accounting  for  QS1 
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This  gives  the  following  integral: 
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To  normalize  we  need  to  perform  the  following  operations 
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Since  we  have  already  normalized  as  a  function  of  angle 
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The  normalized  phase  function  before  integration  can  therefore  be  written  as: 
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The  following  expression  will  be  normalized  to  1  when  integrated  over  angle  and  an  inverse  power  law 
from  0  to  infinity.  It  is  convenient  to  evaluate  related  parameters  such  as  the  asymmetry  parameter  g 
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Note  that  the  first  term  in  the  expression  for  the  cosine  in  terms  of  our  angular  variable  u  is  unity 
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Because  of  the  normalization  of  our  phase  function  expression  we  do  not  need  to  integrate  that  term  as  the 
result  will  be  1.  We  will  now  concentrate  our  efforts  on  the  second  term. 

Let’s  first  perform  the  angular  integral. 

We  obtain: 
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Integrating  over  x  is  analytic  but  leads  to  a  complex  result  involving  2Fj  functions. 
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An  approximate  form  can  be  obtained  directly  from  the  exact  result  by  expanding  the  2Fi  Hypergeometric 
functions  in  series  for  small  arguments. 
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The  formulas  can  be  simplified  by  using  the  following  variable  replacement: 
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And  the  two-term  approximation  becomes: 
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The  rest  of  the  series  terms  can  be  added  if  required  as  follows: 
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This  expression  is  actually  valid  for  the  complete  Fournier-Forand  function  since  the  contribution 
to  the  mean  cosine  of  the  second  additive  term  is  identically  zero  because  it  is  symmetrical 
about#  =  n!  2 . 
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In  order  to  parametrize  the  Foumier-Forand  function  in  terms  of  the  mean  cosine,  we  need  to  reduce  the 
expression  above  to  a  single  parameter.  To  do  this  we  can  use  the  relationship  found  by  Mobley  et  al 
between  index  and  Junge  power  law.  This  relationship  can  be  approximated  as: 

(//  —  3)  =  6  («  —  l) 


Note  that  the  Foumier-Forand  formula  is  only  strictly  valid  between  the  limits  3  <  [l  <  5  ,  which  imply 
an  absolute  limit  in  index  of  0  <  «  —  1  <  1/3 .  However  one  should  be  somewhat  conservative  and  restrict 
the  range  so  as  to  not  approach  the  limits  of  validity  too  closely.  This  issue  should  be  looked  at  in  more 
detail  but  a  suitable  choice  for  the  limits  for  now  would  be  3.25  <  [l  <  4.75 
If  we  substitute  the  Mobley  expression  in  the  equation  for  the  mean  cosine  we  obtain: 
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We  can  easily  obtain  a  simple  approximate  expression  for  (n  —  1 )  in  terms  of  (l  —  g)  by  performing  a 

modified  variable  power  Pade  first  order  approximation  and  inverting  the  resulting  function.  By  this 
method  we  first  obtain  a  formula  valid  from  an  index  of  1.01  to  1.25  which  corresponds  to  a  range  for  g  of 
0.0002  <(l-g)<  0.6: 
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Inverting  this  formula  we  obtain: 
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We  can  now  use  Mobley’s  relation  along  with  this  formula  to  rewrite  the  Foumier-Forand  function  in 
terms  of  the  mean  cosine.  We  simply  need  to  perform  the  following  replacements  in  the  variables  of  the 
FF  function. 
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This  gives  a  form  suitable  for  use  with  the  multiple  scattering  formulas  recently  derived  and  numerically 
verified  by  Piskozub  and  McKee6.  They  show  that  for  order  m  of  multiple  scattering  with  a  single 
scattering  albedo  of  1  the  mean  cosine  is  given  by: 

m 

&  m  =  g 

They  also  show  that  the  mean  cosine  of  the  asymptotic  distribution  for  a  substance  with  a  single  scattering 
albedo  CO  is  given  by: 

_  g(i-ffl) 

We  can  therefore  use  the  parameterized  version  of  the  Foumier-Forand  to  compute  the  multiply  scattered 
distribution  but  we  need  to  keep  in  mind  that  the  limits  on  g  imposed  by  the  range  of  validity  of  the  phase 
function  also  apply  to  both  gm  and  gms . 

The  approach  above  needs  to  be  tested  further  since  there  is  no  evidence  at  this  time  that  the  Foumier- 
Forand  function  is  self  similar  under  the  iterated  finite  convolution  operation  required  to  compute  multiple 
scattering.  As  examples,  the  Gaussian  is  a  self-similar  function  under  convolution  over  an  infinite  range 
and  the  Flenyey-Greenstein  function  has  been  proven  to  be  self-similar  under  the  multiple  scattering 
convolution  operations.  It  is  therefore  highly  unlikely  that  the  final  asymptotic  stationary  state  under  high 
albedo  conditions,  even  though  it  must  have  the  appropriate  gm  and  gms ,  be  described  by  a 

parameterized  FF  function.  The  full  problem  of  going  from  the  parameterized  single  scattering  form  to  the 
correct  asymptotic  distribution  is  still  open  and  will  be  the  subject  of  further  work.  We  now  have  the 
possibility  the  Henyey-Greenstein  function  could  be  used  as  the  asymptotic  multiple  scattering  limit  of  the 
Fournier-Forand  function  with  the  single  scattering  mean  cosine  given  by  the  Foumier-Forand  form 
derived  above.  However,  under  the  moderate  albedo  conditions  encountered  in  ocean  waters,  the  present 
approach  should  give  accurate  results  in  the  important  and  difficult  to  evaluate  transition  range  of  1  to  1 0 
scattering  events. 
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